December 10, 2009 1:42 Applicable Analysis apa-maruno-rcvisc-final 



Applicable Analysis 

Vol. 00, No. 00, July 2009, 1-17 



RESEARCH ARTICLE 

The discrete potential Boussinesq equation and its multisoliton 

solutions 



ON 

o 
o 

(N 
O 

Q 



d 



> 
o 
o 
oo 

od 
O 

o 



Ken-ichi Maruno a * and Kenji Kajiwara b 

^Department of Mathematics, The University of Texas-Pan American, Edinburg, TX 
78539-2999; b Faculty of Mathematics, Kyushu University, 744 Motooka, Fukuoka 

819-0395, Japan 

{December 10, 2009) 

An alternate form of discrete potential Boussinesq equation is proposed and its multisoliton 
solutions are constructed. An ultradiscrcte potential Boussinesq equation is also obtained from 
the discrete potential Boussinesq equation using the ultradiscretization technique. The detail 
of the multisoliton solutions is discussed by using the reduction technique. The lattice potential 
Boussinesq equation derived by Nijhoff et al. is also investigated by using the singularity 
confinement test. The relation between the proposed alternate discrete potential Boussinesq 
equation and the lattice potential Boussinesq equation by Nijhoff et al. is clarified. 
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1. Introduction 

In this article, we propose an alternate form of discrete potential Boussinesq (BSQ) 
equation 



jTm+2jTm—l , 
u n-l u n-l ( 



-hU™~ 1 + 5 2 U™ +2 ) = U^U^i-d.U^ 1 + 5 2 U™\ 2 ) , (I) 



and study the relation to the lattice potential BSQ equation proposed by Nijhoff 
et al. 0, [| 
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Here U = U™ and u = u™ are the dynamical field variables at the site (n, m) of a 
rectangular lattice, and 5±, 5 2 , p, q are the lattice parameters. These are integrable 
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2 The discrete potential Boussinesq equation 

discrete analogues of the potential BSQ equation 

3w tt + 4:C w xx - 6w x w xx - w xxxx = , (3) 
which leads to the BSQ equation^, 0] 

Sutt + 4:C u xx - 2>{u 2 ) xx - u xxxx = , (4) 

by u = w x . Note that Eq.([2]) is the lowest order member of a hierarchy of lattice 
equations which is called the lattice Gelfand-Dikii(GD) hierarchy See also recent 
works about the lattice potential BSQ equation (0)0,0]. 

In this paper, we propose an alternate form of discrete potential BSQ equation 
and present multisoliton solutions. It is shown that the multisoliton solution for 
the discrete potential BSQ equation can be constructed from one for the Hirota- 
Miwa (discrete KP) equation using the reduction technique. Using the proposed 
discrete potential BSQ equation, we can construct the ultradiscrete potential BSQ 
equation. We also study the relation between the alternate discrete potential BSQ 
equation dU) and the lattice potential BSQ equation ([2]). Bilinear equations of the 
lattice potential BSQ equation can be derived systematically using the singularity 
confinement (SC) test. This reveals the relationship with other discrete potential 
BSQ equations. 

2. Multisoliton solutions for the Boussinesq equation 

First, we review fundamental results about reductions of the Kadomtsev- 
Petviashvili (KP) equation. 
It is well known that the solutions of the KP equation 0] 

(-Au t + 6uu x + u xxx ) x + 3u yy = , (5) 

can be expressed via a tau function t(x, y, t) as u(x, y, t) = 2^ log t(x, y, t) where 
r(x,y,t) satisfies Hirota's bilinear equation 

(Dt-AD x D t + 3D 2 y )T-T = 0, (6) 

where D x , Dt, D y are the Hirota D-operators. Soliton solutions of the bilinear equa- 



tion can be written in terms of the Wronskian determinant [8|, |£J, lid . 1 1 11 ] 



T(x,y,t) = Wr (/!,••• J N ) = det(fi n '-%< n>n ,< N , (7) 

with fn^ = f n /dxi , and where /i, . . . , /jv are a se t of linearly independent solu- 
tions of the linear system 

dy dx 2 ' dt dx 3 ' 

For example, ordinary iV-soliton solutions are obtained by taking f n = e® 2 "- 1 +e d2n 
for n = 1, . . . , N, where 

d m {x,y,t) = k m x + k 2 n y + k' i n t + e m -Q, (9) 
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for m = 1, . . . , 2N, where the AN parameters &]_<•••< A^tv and 9i-o, . . . , #2iV;0 are 
real constants. For N = 1, one obtains the single-soliton solution of KP equation: 

Ui,j(x,y,t) = - fci) 2 sech 2 [-(0i - ^)] , (10) 

where i = 1 and j = 2. The most general form of the A^-soliton solution is given 
by 

r(x,y,t) = det(AQK) = ^ Kr> lr .,m„ ^ mi ,.,m» exp mi) ..., mjv , (11) 

l<mi<---<mjv<A^ 

where ^4 = (a n . m ) is the N x M coefficient matrix, = diag(e 01 ,--- ,e 0M ), 
^ m ,,...m = flp, + ••• + Q m t ,, and the M x N matrix if is given by K = 
{k^ 1 ) }^,M,M,M,M,mm^- V mu ..., mN is the Vandermonde determinant 
V mu ...,m N = rii<j<j'<jv(^ - ^) , and A mu „^ m , N is the N x iV-minor whose 
n-th column is respectively given by the m n -th column of the coefficient matrix 
for n = 1, . . . , N. The only time dependence in the tau function comes from the 
exponential phases mi) ..., mjv . Also, for all G € GL(N, IR), the coefficient matrices 
A and A' = G A produce the same solution of the KP equation. Thus without loss 
of generality one can consider A to be in row-reduced echelon form (RREF). One 
can also multiply each column of A by an arbitrary positive constant which can be 
absorbed in the definition of 6q ; o, • • • , #m ; o- 

Real nonsingular (positive) solutions of the KP equation are obtained if k\ < 
■ ■ ■ < kM and all minors of A are nonnegative. Under these assumptions and some 
fairly general irreducibility conditions on the coefficient matrix, Eq. (|23p produces 
(N-, A r + )-soliton solutions of the KP equation with N- = M — N and N + = N, as 
in the simpler case of fully resonant solutions. Asymptotic line solitons are given 
by Eq. (ll0p with the indices i and j labeling the phases Oi and 8j being swapped 
in the transition between two dominant phase combinations along the line Oi = Oj. 
Asymptotic solitons can thus be uniquely characterized by an index pair with 
1 < * < j < M. Recently, line soliton solutions of the KPII equation were clas- 
sified using this formulation [l^, EH, [3, EEI, E3, G3] • Elastic 2-soliton solutions are 
classified into three classes: ordinary (O-type), asymmetric (P-type) and resonant 
(T-type). The coefficient matrices corresponding to these classes have the following 
RREFs: 

. Aa00\ , /l00-6\ , /10 -c-d\ , 1oN 
A °={oOlb)> Ap =[01a } ^ T= (oi a b > < 12 > 

where a,b,c,d > are free parameters with ad — be > 0. These three types of 
solutions cover disjoint sectors of the 2-soliton parameter space of amplitudes and 
directions. Moreover, their interaction properties are also different. This difference 
is obvious in the case of T-type solutions, but also applies to O-type and P-type 
solutions, since P-type solutions only exist for unequal amplitude, and the interac- 
tion phase shift has the opposite sign for O-type and P-type solutions. We remark 
that inelastic 2-soliton solutions fall into four categories, but we do not focus on 
them in this article. 

It is known that the KP equation reduces to the Korteweg-de Vries (KdV) equa- 
tion 



- Au t + 6uu x + u xxx = 



(13) 
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by the symmetry constraint du/dy = 13, l2Q . l2ll . [22j. In the bilinear form, the 
bilinear KP equation ([6]) is reduced to the bilinear KdV equation (D x — 4D x Dt)r ■ 
t = by the constraint of omitting terms including D y , which is the so-called 2- 
reduction. As mentioned in [l^, this constraint implies that all the solitons of the 
KdV equation are parallel to the y-axis. Thus we must have a condition kj = — ki 
for each [i, j']-soliton. From the ordering k\ < hi < ■ ■ • < &2JV, we must assume 



k\ < k2 < • • • < fcjy < , kjy + j = —kjy-j+i for j = 1, 
This allows only P-type soliton solution in which ^4-matrix is 



,N. 



A 



/l ■■■ 
1 ••• 

\0 • • • 1 aj\r,jv • • • 



ai,27v\ 

0>2,2N-1 







J 



Let us consider multisoliton solutions of the BSQ equation. It is also known that 
the KP equation is reduced to the BSQ equation without the term u xx , namely 



35 u t 't' ~ 3(/u 2 ) 



xx ^^xxxx 



o, 



(14) 



20, 21, 2 



by the symmetry constraint du/dt 
independent variable t' such that y = y/—5ot', So 
bilinear KP equation ([6]) is reduced to (D x — 35oD 2 )t 
of omitting terms including D t , which is the so-called 3-reduction. This condition 
implies that all the solitons of Eq. (ll4|) are parallel to the i-axis. Thus we must have 



. Here we introduced a new 
±1. In the bilinear form, the 
t = by the constraint 



a condition k 2 + kjki + kf = 0, i.e. kj = ujki or kj = u! 2 ki (uj = —1/2 + i\/3/2, 
u> 2 = —1/2 — i\/3/2, a; 3 = 1), for [i, j]-soliton. In the BSQ equation, k\,...,k2N 
take complex values in general. Since k\,...,k2N are complex values, we cannot 
consider the ordering of ki, A^jv which was assumed when we consider KP soliton 
solutions. However, from the constraint of kj, we have a restriction to the ^-matrix 
such that each row has only 2 or 3 nonzero elements and each column has only 
one element. Note that there are only 2 nonzero elements in each row and only one 
element in each column in the case of the KdV equation because of the constraint 
kj = —ki 131 ] . Note also that there is no such restriction in the case of the KP 



equation. For 2-soliton solutions, 2 types of 2-soliton interactions are possible, i.e. 
2 elastic soliton interactions (O-type and P-type), other 2-soliton interactions are 
impossible because some columns in A-matrix have 2 elements. However, there is 
no distinction between O-type and P-type solitons because we do not have the 
ordering of k±, ...,k2N- Thus 2-soliton interaction of the BSQ equation is actually 
only of one type. To get real solutions, we must remove imaginary numbers by the 
following ways. 

The case of 2 elements: 
(i) Suppose that we have 2 elements in the i-th row of the ^4-matrix. Let the 
corresponding wave numbers of these elements be kj 1 , kj 2 ; (ii) Let kj 2 = ukj 1 
(or kj 2 = u} 2 kj t )] (hi) Using the gauge invariance of r-function, each element in 
the Wronskian determinant can be fi ~ 1 + ajj 2 exp((c 
l)k 2 t' + 6j 2 fl) (or fi ~ 1 + aij 2 exp((uj 2 - l)kj t x + VSq(u> 
(iv) Reparametrize kj 1 and kj 2 by Kj = (uj — l)kj 1 and = y— oo 
Ki = (u) 2 — l)kj 1 and = \J—5q(uj — 1)^). Then we have fi ~ 1 + a 
Qit' + 6j 2 fl) with linear dispersion relations k\ — 3S^0\ = 0; (vi) So a set of /, 



l)k j± x + y^o^ 2 - 
l)k 2 t' + 9 nfl )); 

M" 2 - l)kl (or 

i,j 2 exp(/tjX + 
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1 + dij 2 exp(KjX + Ojt' + 9j 2 fl) with linear dispersion relations nj — 35oQ 2 = gives 
the real and regular multisoliton solutions. 

The case of 3 elements: 
(i) Suppose that we have 3 elements in the i-th row of the ^4-matrix. Let the 
corresponding wave numbers of these elements be kj ± , kj 2 , kj 3 ; (ii) Let kj 2 = u>kj 1 
and kj 3 = ui 2 kj 1 ; (iii) Using the gauge invariance of r-function, each element in the 
Wronskian determinant can be fi ~ l + aj j2 exp((a; — l)kj t x + (ui 2 — l)k 2 i y + 9j 2 fi) + 
dij 3 eyLp((tL> 2 — l)kj 1 x+(uj—l)k'j i y+6j 3 fl); (iv) Reparametrize parameters kj 1 , kj 2 and 
kj 3 by K it 2 = kj 2 -k h , K ij3 = k j3 -k jl , Q i>2 = x/^^fc? -fc?) and Q i)3 = \^S^(k 2 3 - 
k 2 J. Then we have fi ~ 1 + aij 2 exp^x + fiji 7 + #j 2i0 ) + ajj 3 exp(KjX + il.it' + 9j 3: o) 
with linear dispersion relations nfj — 35o^l 2 j = for j = 1, 2; For example, consider 
N = 1. This gives Y-shape soliton resonance interaction. 

The KP equation is reduced to the BSQ equation 

35 u t >t' + 4c u xx - 3(u 2 ) xx - u xxxx = , (15) 

by the symmetry constraint = cojj- Here we introduced a new independent 
variable t' such that y = y/—5ot'. In the bilinear form, the bilinear KP equation 
([6]) is reduced to (D x — 4cqD 2 — 35qD 2 )t ■ r = by the constraint of replacing D t 
by cqD x , which is the so-called 3-pseudo reduction [12]. To realize this constraint 
in the multisoliton solutions, we can add the constraint kj — kf = co(kj — ki), i.e. 

kj = o ( — ki ± \ / 4cq — 3kf ) , to the soliton solutions of the KP equation. Note that 



*j — 2 

kj can be real if 4co > 3k 2 . So we can assume the ordering of k\ < .. . < k 2 ^ ■ From 
the constraint we have a restriction to the A-matrix such that each row has only 2 or 
3 nonzero elements and each column has only one element. For 2-soliton solutions, 
2 types of 2-soliton interactions are possible, i.e. 2 elastic soliton interactions (O- 
type and P-type), other 2-soliton interactions are impossible because some columns 
in ^4-matrix have 2 elements. 

Let us consider explicit real and regular multisoliton solutions. For the case of 
Jo = —1, (i) Suppose that we have 2 elements in the i-th. row of the ^4-matrix. 
Let the corresponding wave numbers of these elements be kj t , kj 2 ; (ii) Let kj 2 = 

\ {^kji i y^4co — ; (iii) Each element in the Wronskian determinant can be 

fi ~ a i,jx exp(kj 1 x+k 2 i t'+9j 1 fi)+aij 2 exp(kj 2 x+k 2 2 t' +9j 2i q). This gives multisoliton 
solutions for the BSQ equation. 

If we have 3 elements in the i-th. row of the ^4-matrix, each element in the 
Wronskian determinant can be fi ~ ex.p(kj 1 x + k 2 t' + 9j l: o) + aij 2 ex.p(kj 2 x + 
k 2 2 t' + 9j 2 fi) + ciij 3 exp(kj 3 x + k 2 t' + 9j 3t o). For N = 1, this gives Y-shape resonant 
soliton solution. 

For the case of 5$ = 1, (i) Suppose that we have 2 elements in the i-th row of the 
^-matrix. Let the corresponding wave numbers of these elements be kj 1 , kj 2 ; (ii) Let 

kj 2 = i y— kj 1 ± yj 4co — 3/c^ ^ ; (iii) Using the gauge invariance of r-function, each 
element in the Wronskian determinant can be fi ~ l + aj j2 ex.p((kj 2 — kj^x + '^k 2 ^ — 
kjjt' + 9j 2 fl); (iv) Reparametrize kj 1 and kj 2 by Ki = kj 2 —kj 1 and $7^ = i(/c| 2 — kjj. 
Then we have fi ~ 1 + a,ij 2 exp^rr + fijt' + 9j 2 $) with linear dispersion relations 
k\ - Acqk 2 - 35 Q 2 = 0. 

Next, we consider the case of complex parameters. Suppose that a parameter ki 
corresponds to the pivot of A-matrix. If 4cq < 3kf, then we have a complex param- 



eter kj = 2 y— h ± iy 3kf — 4coJ among k\, k 2 N ■ Since some of k\, k 2 N are 
complex values, we cannot consider the ordering of fci, k 2 N which was assumed 
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when we consider KP soliton solutions. However, from the constraint we have a 
restriction to the ^4-matrix such that each row has only 2 or 3 nonzero elements 
and each column has only one element. For 2-soliton solutions, 2 types of 2-soliton 
interactions are possible, i.e. 2 elastic soliton interactions (O-type and P-type), 
other 2-soliton interactions are impossible because some columns in ^4-matrix have 
2 elements. However, there is no distinction between O-type and P-type solitons 
because we don't have the ordering of k\, feiv- Thus 2-soliton interaction of 
the BSQ equation is actually only one type again. To get real solutions, we must 
remove imaginary numbers by the following ways. 

The case of 2 elements: 
(i) Suppose that we have 2 elements in the i-th row of the ^-matrix. Let 
the corresponding wave numbers of these elements be kj lt kj 2 ; (ii) Let kj 2 = 

\ ( — kji i iy/3/c^ — 4co^ ; (iii) Using the gauge invariance of r-function, each el- 
ement in the Wronskian determinant can be ~ 1 + ajj 2 exp((kj 2 — kj 1 )x + (fe| 2 — 
hji)y~^~@h,o)> Reparametrize kj t and kj 2 by Ki = kj 2 — kj 1 and = \/— 5o(fc| 2 — 
k^). Then we have fi ~ 1 + Oj j2 expire + flit' + 9j 2t o) with linear dispersion rela- 
tions k\ — AcQnf — 35o^f = 0; (vi) So a set of fi ~ 1 + Oj j2 exp(KjX + Q{t' + 0j 2 $) 
with linear dispersion relations k\ — 4cok| — 30? = gives the real and regular 
multisoliton solutions. 

The case of 3 elements: 
(i) Suppose that we have 3 elements in the i-th row of the ^4-matrix. Let the 
corresponding wave numbers of these elements be kj 1 , kj 2 , kj 3 ; (ii) Let kj 2 = 

\ (~% + l \f^ji ~ ^ c °) anc * % = \ ~~ iy/3/c^ — 4co^ ; (hi) Using the gauge 

invariance of r-function, each element in the Wronskian determinant can be 
fi ~ 1 + a ij2 exp((fc j2 - kj^x + v /Z ^o(^ 2 ~ tfjt' + %, ) + a ij3 exp((fcj 3 - kjjx + 
V— 5o(A;| 3 — kjjt' + Oj 3 fl); (iv) Reparametrize kj t , kj 2 and /cj 3 by = % 2 — kj t , 
«i,3 = % ~ %u fi i,2 = v / -^o(^] 2 - fcfj and O i;3 = </^(/c| 3 - k?J. Then we 
have fi ~ 1 + ajj- 2 exp(KjX + fiji' + 0j 2 ,o) + Ojj 3 expire + Oji' + 9j 3t o) with linear 
dispersion relations re^ ■ — 4coftfj — 3<5o^fj = for j = 1, 2; For example, consider 
N = 1. This gives Y-shape soliton resonance interaction. 



3. Discrete analogues of the potential Boussinesq equation 

3.1. The discrete potential Boussinesq equation 
Here, we present the main theorem in this article. 
Theorem 3.1 : The difference- difference equation 

uZFuZJi-Sivr 1 + <W +2 ) = ^ + X m -iMi£Ci + hu™\ 2 ) , (16) 

where Si = 02(01 — 03) and 82 = 03(01 — 02) and a±, 02, 03 are arbitrary real 
constants, is an integrable discrete analogue of the potential BSQ equation |3[). 
Moreover, the discrete potential BSQ equation has multisoliton solutions 

T m 

U n n = J ^, T™ = det(ASP), (17) 

'n 

where A = (ot n%m ) is the N x 2N coefficient matrix, 6 = diag(e 6 * 1 , • • • , e e2N ), e * = 
Pj(l—pjai)~ n (l—pja2)~ m , and the 2NxN matrix P is given by P = (p 1 ^ 1 )- where 
the 2N parameters pi,...,P2N are rea ^ constants. The A-matrix has a restriction 
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such that each row has only 2 or 3 nonzero elements and each column has only 
one element. In the case having 2 elements and (i,j 2 ) in the i-th row of the 

A-matrix, pj 2 must satisfy a reduction condition 



J_ (1 - azPfr) ygglj ~ Q3Pji)(Q2 ~ 4a 3 + 3a 2 a3Pj 1 ) , . 

Pj2 ~ ' n n ' \ °) 

a 2 2a 3 2,a 2 a 3 

Proof : The Hirota-Miwa (discrete KP) equation is written as 

oi(o 2 - a 3 )r(ni + 1, n 2 , n 3 )r(m, ra 2 + 1, ra 3 + 1) 
+a 2 (a 3 - ai)r(ni, n 2 + 1, n 3 )r(ni + 1, ra 2 , n 3 + 1) 
+a 3 (ai - a 2 )r(n 1 ,n 2 ,n 3 + l)r(ni + 1, ra 2 + 1, n 3 ) = , (19) 

where r depends on three discrete independent variables n\, n 2 and n 3 , and a\, a 2 
and a 3 are the difference intervals for n\,n 2 and n 3 , respectively [13]. 
The Casorati determinant solution for the Hirota-Miwa equation (|19|) is as 



follows [24| : 



r(ni,n 2 ,n 3 ) = det(^(?ii, n 2 , ra 3 ; s + j - l))i<ij<jv , (20) 
where ipi, ipw are a set of linearly independent solutions of the linear system 

A n .ipi(ni,n2,n 3 ;s) = ^(ni,n 2 ,n 3 ; s + 1) , (j = 1,2,3). 
Here A nj are the backward difference operators: 

X.JU,) f{ "' ' j . (,- = 1,2,3). (21) 

For example, ordinary iV-soliton solutions are obtained by taking 

ip i (n 1 ,n 2 ,n 3 ;s) = a 2 i-ip 2 i-i0- -P2i-iai) _7lI (l - p 2 j_ia 2 )~ n2 (l - f> 2i -ia 3 )~ n3 
+a 2 ^(l -p 24 ai)- ni (l -p 2i a 2 )- n2 (l -p 2i a 3 )- n3 , (22) 

for n = 1, ...,iV where the 4N parameters p\ < ■■■ < p 2 jy and a±,--- ,a 2 N are 
positive real constants. The most general form of the iV-soliton solution is given 
by 

r(ni, n 2 , n 3 ) = det(„4GP) = ^ Kn 1; ...,m w Am,...,m w ex P , miv , 

l<mi<--<m JV <2A r 

(23) 

where A = (a n ,m) is the N x 2N coefficient matrix, = diag(e ei , • • • ,e d2N ), 
e d s = p*(l - pja!)-™^! - p ia2 )~ n2 (l - pja 3 )- n3 , and the 2iV x N matrix P is 
given by P = (p™ -1 ). V^. lr .. )mjv is the Vandermonde determinant V mu _ >mN = 
rii<3<j'<Jv(P'" 3 -; — Pm 3 ) i and ^4 mi) ..., mjv is the TV x iV-minor whose n-th column is 
respectively given by the m n -th column of the coefficient matrix for n = 1, . . . , N. 
For all G G GL(N, IR), the coefficient matrices A and .4/ = G A produce the same 
solution of the Hirota-Miwa equation. Thus without loss of generality one can 
consider A to be in RREF. 
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Let us consider the 3-reduction condition 

r(ni, n 2 , n 3 ) o r(m, n 2 - 2, n 3 - 1) , 



(24) 



for arbitrary ni,n 2 and re 3 Applying this reduction condition, we can omit the 
dependence on n 3 and obtain the bilinear form 

ai(a 2 - a 3 )r(ni + 1, n 2 )r(ni, n 2 - 1) + a 2 (a 3 - ai)r(ni,n 2 + l)r(ni + l,n 2 - 2) 
+a 3 (ai - a 2 )r(ni, n 2 - 2)r(ni + 1, n 2 + 1) = . (25) 

After the change of variables n\ — > n, n 2 — > m, r(ni,ra 2 ) — > r™, we obtain 

ai(oa - a^V™ + a 2 (a 3 - aOC+V^ 1 + a 3 ( ai - a^T™" 1 ^ 2 = , 

(26) 

which is the bilinear form of the discrete potential BSQ equation. 

Now we impose a constraint on the parameters of the solution so that the reduc- 
tion condition is satisfied. For simplicity, we consider the case in which ipi, 
have 2 terms. Then we observe 

ipi{ni,Ti2 + 2,113 + 1; s) 

= p s jl (l-p jl a 1 )-^(l-p jl a 2 )-^- 2 (l-p jl a s )- n -- 1 
+p%(l- Pj2 a 1 )-^(l-p j2 a 2 y^- 2 (l-p j2 a 3 )-^- 1 

= pl(l-p jl a 1 )- n Hl-p jl a 2 )- n *- 2 (l-p jl a 3 T n3 - 1 



1 ~Pj2 a l 

1 -p 31 ai 



where 



1 -p j2 a 2 



1 -p n a 2 

If we apply the reduction condition 



1 — pj 2 a 2 
l-p n a 2 



-2 /, x -1 

1- Pjl a 3 



l-p j2 a 3 
l-p n a 3 



(1 -p h a 2 f{l -Pj^az) = (1 -p h a 2 ) 2 (l - p j2 a 3 ) 



(27) 



(28) 



i.e. 



1 (1 - qggjj 1/02(1 - a 3 Pj 1 )(a 2 - 4a 3 + 3a 2 a 3 p jl ] 



a 2 2a 3 2a 2 a 3 

we obtain 

ipi(ni,n 2 ,n 3 ; s) = (1 - p, 2 a 2 )~ 2 (l - Pj 2 a 3 )~ 1 ip i {n 1 , n 2 - 2,n 3 - 1; s) 
Finally we have a reduction condition 

iV 

r(m,n 2 ,n 3 ) = JJ(1 -Pka 2 )~ 2 (l - Pka 3 )~ 1 T(n ll n 2 - 2,n 3 - 1) . 
fe=i 



(29) 



(30) 



(31) 



December 10, 2009 1:42 Applicable Analysis apa-maruno-rcvisc-final 



Applicable Analysis 9 

In a similar way discussed in section 2, we can consider the general multisoliton 
solutions of the discrete potential BSQ equation. □ 

The reality condition is pj 1 < 1/03 , pj 1 > (—02 + 4a3)/(3a2a3) for 02 < 03, and 
Pj 1 > 1/03, pj 1 < (— 02 + 4a3)/(3a2a3) for 02 > 03. With parameters satisfying 
these conditions, we can construct real and regular multisoliton solutions using the 
formula in Theorem 13.11 

Otherwise, we must use the technique of reparametrization which was used in 
section 2. The procedure is as follows, (i) Suppose that we have 2 elements in 
the z-th row of the A-matrix. Let the corresponding complex wave numbers of 
these elements be pj 1 and pj 2 ; (iii) Using the gauge invariance of r-function, each 
element in the Wronskian determinant can be fi ~ 1 + aij 2 K^ n £l^ m where Ki = 
(1 — a\Pj 2 )/(l — aipjj, Oj = (1 — a2Pj 2 )/(l — a,2Pj 1 ); (iv) Choose new parameters 
Ki and Oj to be real numbers. 

It is easy to take the ultradiscrete limit in Eq. (|16p . 

Theorem 3.2 : The ultradiscrete potential BSQ equation is 

V™? + V^ 1 + L™+ 2 = V™ +1 + C-i + K¥ - (32) 
U n m+2 = max(L™+ 2 , K"" 1 + ci ) - c 2 . (33) 

Proof : Use the standard procedure of ultradiscretization. Equation is rewrit- 
ten in the form of 



TTm+2jjm—l jm+2 



'n-l 
rm+2 



TTm+ljTm rm+2 
u n {J n-\ 1 n~\ > 



-hU™- 1 + 5 2 UT +2 



(34) 
(35) 



Introduce new variables U™ = exp(V^ 7l /e), I™ = exp(L™/e), 5\ = exp(ci/e), 82 = 
exp(c2/e). Then take the limit e — > + using the formula lim e ^o+ eln(exp(A/e) + 
exp(B/e)) = max(A, B) for A, B 6 R>ol3- □ 

Remark 1: 

Applying the gauge transformation U™ — ► (—l) n+m U™ to Eq. (fT6]) . we have 



jTm+2jjm—l rm+2 



Tjm+ljjm, rm+2 
u n ^n-l-'n-l ' 



rm+2 



m+2 



(36) 
(37) 



This form gives a more symmetric form of the ultradiscrete BSQ equation 



T/m+2 . i/m-1 1 rm+2 _ T/m+1 , r/m , rm+2 /oo\ 
V n-1 + V n-1 + L n ~ V n + V n-1 + L n-1 ' \ 6b > 

L™ +2 = max(U n m - 1 + ci , U n m+2 + c 2 ) • (39) 



Remark 2: 



Date et al. proposed another discrete potential BSQ equation[26jj 



,m— 1 



J n-l u n 



(ai(o 2 - a 3 )v™ +1 + a 2 (a 3 - ai)u™ +1 ) 



C-i^n-tfM^ - as)^" 1 + «2(«3 - aOd; 



(40) 
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By the transformation v™ = t™_^~i /t™, we obtain a bilinear equation 

ai (a 2 - a 3 )C +1 r™ ! + a 2 (a 3 - a^r^C -1 + a 3 (ai - a 2 )r^r^ = . 

(41) 

This bilinear equation is obtained by adding the reduction condition 

r(rei, n 2 , n 3 ) o r(m - 1, n 2 - 1, n 3 - 1) , (42) 

which gives yet another 3-reduction. 

For this discrete potential BSQ equation, we can also make an ultradiscrete 
analogue of the potential BSQ equation 



A n-1 + A n + r n ~ A n-1 + A n+1 + Y n-1 i l 4,: V 



L n+1 

„, _„/ x^rn— 1 



= max(X-Y + C1 ,X- + 02) , (44) 



taking the ultradiscrete limit after setting v™ = exp(X™/e), w™ = exp(Y™/e), 
a\ = exp(ci/e), a 2 = exp(c 2 /e) where w™ = aiv™^ + a 2 v™, ax = a\(a 2 - a 3 ), 
a 2 = a 2 (a 3 - ax). 



3.2. TTie lattice potential Boussinesq equation 

Singularity confinement (SC) test was proposed by Grammaticos et al. as a detector 
of integrability in discrete systems 271 ] . This property has been applied to several 
problems [H, [2^, [13, HH, HI]. The SC test is also powerful tool for constructing 
solutions for discrete integrable systems (30l. 31, H2]- In this section, we apply the 



SC test to the lattice potential BSQ equation of Nijhoff et al. and obtain bilinear 
equations using the result of SC test. 
We start from the slightly simplified form of the lattice potential BSQ equation 



P 



P 



1 + a n+X 



U' 



{p + 2q + u™ +1 
-(p + 2q + u™ 



P~ 

,m+2 
l n+2 



q + < +2 



u, 



m+1 
n+X 



q + u™+ 2 - uT + 



n+X 



)(p-q + u, 



n+X 
m+l 



*n+2 



U 



n+X; 



(45) 



Introducing new variables 



u, 



V 



m-2 



p + 2q + u™Ii 



(46) 



Eq.(]45j) is written as 

„3 



P ~ Q 
1 n+1 



p" - q" 



jin+X 
1 n+X 



-irm+2 jm+X 
V n+2 1 n+2 



-irm+2 jm 
V n+X J n+1 ' 



jm I T/m+2 _ jm+z , -trr 
1 n V n+X — 1 n+X v n 



m+2 



rrn+3 



(47) 
(48) 



After the transformation of the independent variables m + n 



m, Eas.diSD.diBI) 
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and (|47h . (I48h are rewritten in the following form: 

p3 _ g3 p3 _ g3 



q + <ti - <ti p-i + < +2 - <+ 2 

= (p + 2q + <!+/ " <V + 2 4 )(P " Q + <C + + i 3 " <$) 

-(p + 2q + u™ - <!+ 3 )(p - 9 + C +1 - C+i 1 ) , (49) 



C=P-9 + C_i-C K m =P + 2g + C-"i 3 -<, (50) 



and 



3 3 3 3 

P ~ Q _ P ~ Q _ T/ m+4 T m+3 _ T/ m+3 T m+1 /r-i \ 

jm+2 jm+2 ~ V n+2 1 n+2 v n+l 1 n+l > V oi J 

J n+2 i n+l 

fffl i T/m+3 _ rm+3 , T/m+3 (KO\ 

J n >n+l — 1 n+l *n ■ l oz J 

After some calculation, we get the form which is suitable to perform the SC test: 

V m_ 1 ( ym-l rm-3 P 3 ~ ^ p 3 ~ q 3 \ , , 

J " \ L n 1 n-l J 

C = C-i 3 + K" " V£ x . (54) 

Now we perform the SC test. Suppose that we have some initial data which 
can be evolved using the above system. Let us assume that during the successive 
applications of the above system the value of / at (m — l,n) becomes zero. The 
point where this occurs depends on the initial data. Then we have the following 
pattern of and oo: 

C C+i, Cff > - °°. °} . {V™, C-t 3 } - {oo, oo} . 

One can see nonzero finite values for all dependent variables in further steps, so 
the singularity is perfectly confined. Suppose that this pattern was created by a 
function F™ which has a zero at (m, n). 

Using the above singularity pattern, we obtain the independent variable trans- 
formation 

rra+l Tpm—1 

C = « n FmF ^ ■ (55) 

Since u m and I™ are related by Eq. ([50j) . u m should have F™ in the denominator. 
Thus we set 

rim 

< = ^ ■ (56) 



From Eqs. (j46j) and (|56j) . we obtain 

/im rim— 3 (~i m 

In=P-9 + Py n r 1 -P F ^, V™= P + 2 q + P-^±-p±L. (57) 

n—1 n n—1 n 
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Using Eqs. (f53~l) . (|54|) . ([55]) and (|57|) . we obtain the following equations: 

(p - q)F^_- 1 1 Fr 1 + PG^F™- 1 - PF^G™- 1 = aF™F™~* , (58) 
l—^F™ +1 F™ +l - o(p + 2g)F n m + + 1 2 F n m - 1 

-a/^r 1 ^ 2 + ^r'C = Tic^r+V . ( 59 ) 

pmpm-2 3 _ 3 r-m+1™ ^/F m F m+1 

j n n— 1 _ P </ n— 1 n 1 ' n— 1 n 

™ 1 ™ 1 O ™ I o ™ 1 I 



pm-lpm-1 a 2 F rn+2 F m-l pm-1 

n n—1 n n—1 n n—1 

pm+3pm+l 3 3 pm+lpm pmpm+1 

r n+l r n _ P — Q r n r n+1 7 r n r n+1 / fif) x 

a rpm+2rpm,+2 a 2 pm+2 nra-1 a nra+2 pm-1 ' V DU / 

n+1 n n+1 n n+1 n 



where 7 is a decoupling constant. Note that V is written in the following form: 

rmi-2 (~im+\ 3 _ „3 Z?m pm-1 pm-1 p- 

ym_ ( 1 2fl 1 I rc- 1 " - £ * J 7 » fRi 1 ) 

"n T" T p m _ 2 M F m+l 2 p m+l p m-2 pm+1 pm-2 • 

n—1 n n "n—1 n n—1 

Assuming 7 = 0, we obtain 



(p - q)F n ^ 1 1 F n n - 1 + PG^F™- 1 - 0F™- x l G™- 1 = aF^F^ 2 , (62) 

l^F™ +1 F™ +1 - a(p + 2g)F™+ 2 iC- 1 - 1 f£+ 2 

+a/3F n m - 1 G™+ 2 = 0, (63) 
3 3 

pm-lpm+3 P ~1 pm+2pm _ zpmpm+2 /r>A\ 

"^n r n+l a 2 n r n+l ~~ " r ™ r n+l ' v U / 

where (5 is a decoupling constant. After changing back to original independent 
variables (m — n — > m) , we have 

(p - q)F^_ 1 F^ + PG^F™- 1 ~ W-iCr 1 = ai^l^ 1 , (65) 

l^F^F^ 1 ~ «(P + 2g)F™V iC" 1 - a^G^F^ 

+a/3F n m - 1 G-+ 1 = 0, (66) 

3 3 

™ jpm-l Tpm+2 Pj~q 77im+2 pm- 1 r 771m pm+1 /(?7\ 

^n+1 ^2 ^" ^n+1 ~ <^n ^n+l ' l 57 ) 



Note that Eq. (|67p can be derived by eliminating G from Eqs. (|65p and (J66J). This 
is a discrete analogue of a bilinear form of the potential BSQ equation (J3j) - Thus 
we have the following theorem. 

Theorem 3.3: Solutions of Eg. (4°V, i.e. Eqs. and $J3$ , are expressed in the 
following form: 

(~rm pm+l pm-1 svm+l ^ m 

m _ jm _ _ n 1 n-1 _„_„(_ n-1 _ 

r n r n r n-l * n-1 n 

rmi-2 f-im 

V m - n -4- la A - — 

r n-l n 
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where F™ and G™ satisfy Eqs. H)5\). [66\) and [67\ ). Moreover, F™ is given by t™ 
of the discrete potential BSQ equation p6\) in Theorem \3.1\ when parameters have 
relations 



a o 3 (oi - a 2 ) p 



3 „3 



a 2 [ai - a 3 ) 



5 ai(a 3 -a 2 )' a 2 5 ax(a 3 — a 2 ) 



Furthermore, we have the following theorem. 



Theorem 3.4 : The above r-functions F™ and G™ are given by 



rpm m /~im m , 

^ n 'n i Pn ' 



13 aij \{3 a 2 



(68) 



t™ is defined by Eg. 169\) (t™ is also given in Theorem 3.1), p™ is defined by 
Eg. (f70| j. Here parameters should be chosen as a 3 = —a 2 /2, 



ai a 2 J a 2 



(p 3 - g 3 )afa| 



(2ai + a 2 )(ai - a 2 ) 2 



Proof : Let ns intro duce the following r-functions using the Freeman-Nimmo no- 
tation |a 012a 1331: 



Pn(k) 



ipi(n,m, k; s) ipi(n,m,k;s + l) ■■■ ipi(n, m, k; s + N — 1) 
ip2(n,m,k;s) ip2(n,m,k;s + l) ■■■ ip2(n, m, k; s + N — 1) 

4>N(n, m, k; s) ipN(n, m,k;s + 1) ■ ■ ■ ipN(n, m, k; s + N — 1) 

0,1,.-- ,N-2,N-1|, (69) 

ipi(n,m,k;s) ■■• ipi{n, m, k; s + N — 2) ipi{n,m,k;s + N) 
ip2(n, m, k; s) ■■■ ip2 (n>, m, k; s + N — 2) tp2 {ti, m, k; s + N) 

ipKi{n, m,k\s) ■ ■ ■ %j)jsi(n, m, k; s + N — 2) ipN(n, m i k; s + N) 
0,1,..- ,N-2,N|, (70) 



where ipi, ipjy are a set of linearly independent solutions of the linear system 



A n ipi(n, m, k; s) = A m ipi(n, m, k; s) = A fc Vi(n, m, fc; s) = ^(n, m,k;s + 1) 



Here A re , A m , Aj, are the backward difference operators 



AJW = '<">- '<-■> , Am/(m)5 /(m)-/(m-l) i A ^/W -/(*-« 



"1 



«2 



0.3 
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Then we notice the following difference formulas [24], |3 



ai-Ci(fc) = |0,1, •■■ ,N-3,N-2,N-2 n+ i|, 
a 2 T™ +1 (k) = ]0,1, ,N-3,N-2,N-2 m+1 |, 
a 3 T™(k + l) = |0, 1, - - - ,N-3,N-2,N-2 k+1 |, 
(ai-a 2 )T^ 1 1 (k) = |0,1, ,N-3,N-2 m+1 ,N-2 n+1 |, 
(fl2 - a 3 )r™ +l (k + 1) = |0, 1, • • • , N - 3, N - 2 fc+1 , N - 2 m+1 | 
aiPn+i(k) ~ r™ +1 (k) = |0, 1, • ■ ■ , N - 3, N - 1, N - 2 n+1 | , 
a 2 p™ +1 (k)-T™ +1 (k) = \0,l,--- ,N-3,N-l,N-2 m+1 |, 
o 3 p™(fc + l)-T™(fc + l) = ]0,1, ,N-3,N-l,N-2 k+1 



From Plucker relations and these difference formulas, we obtain the following three 
bilinear equations 



prwc + iW-^ + i(fc)c +1 w 

= (---) - , (7i) 

\«1 0,2/ 

P ™(k + i)C +1 (fc) - p™ +1 {k)<{k + 1) 

= ( - - -) «+^ k + - + 1)) . (72) 

V°l a 3/ 

ai(a 2 - a 3 )C +1 (A:)C +1 (A : + 1) + a 2 (a 3 - ai )r™ +1 (k)r™ +1 (k + 1) 

+a 3 (a 1 -a 2 )C(fc + l)C+ 1 (fc) =0. (73) 



Applying the 3-reduction condition r™{k + 1)ot™ 2 (/c) and setting f™ = T™(k), 
g™ = p™(k), we obtain 



L im m ir 

Jn+l Un+lJn 



m— 2 fm 
Un Jn+l 



m fm—2 
Un+lJn 



1 



1 

«2 



a m+1 j^m fm+l fm \ 
n+l Jn in Jn+l) 



(74) 



| / ern—2 fm fm fm—2\ (7K\ 
I Un+1 Jn ~~ Jn+lJn ) i I ' °J 



ai(o2 - a 3 )/ n m +1 / n m - i + o 2 (a 3 - a 1 )C +1 f^~ l 2 + a 3 (ai - a,)/™"^ 1 = 



cm— 2 fm+l 



(76) 



To compare these bilinear equations with bilinear equations (|65p - (j67p . we set 



03 



^2 rm 
" 2 ' n 



pn + gm 
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Then bilinear equations (|74j) -(|76 p and (165p -(|67p are rewritten as 

im+l fin fm+1 / - | fm+1 fin {77\ 

n n Jn+1 ViJn ~~ I ~ ~ I Jn+1 Jn > \'< ) 

\«1 «2/ 

? m— 2 fin t m fin— 2 [ . ^ \ rra-2 j-m /'7S\ 

n n Jn+1 ~~ n n+lJn ~ \ I I Jn+1 in i l'°J 

\Ol 02/ 

SaiC+iC" 1 - (2ai + a 2 )/ n m+1 /r+l 2 " («i " a 2 )f™~ 2 f™+ 1 1 = , (79) 

and 

W +1 iC+! " W+l^n m+1 = Offfi 1 ^ , (80) 

o rrm-2 77m Rum rpm— 2 P 9 77m— 2 77m /oi\ 

P-"n ^n+1 - P H n+l*n = ^2 ^"+1 ^ > l 81 J 

Assuming F™ = f™,H™ = h™ and comparing coefficients of these equations, we 
have the relations of parameters 



ax a 2 J a 2 



(j> 6 - q 6 )a{a2 
(2ai + a 2 ){ai - a 2 ) 2 



□ 



Note that Hietarinta and Zhang gave the Casorati determinant form of the 
r-functions F™ and We also note that if we introduce the auxiliary 

independent variable x by d x ipi(n,m,k;s) = ipi(n, m, k; s + 1), then p™ can be 
expressed as p™ = d x T™. 

Remark 3: 

In the case of p 3 — q 3 = 0, Eq. (f59j) is 



- a(p + 2q)F™+ l l F™- 1 - a^- l F^ + af3F^~ 1 G^ = 7*^1 . (83) 
Equation (fBTJl) is written as 



pro 77771- 1 pm+ip+1 pm+2 rpm+1 77m rpm 
r n r n-l . 7 1 n-1 1 n r n+1 r n . 7 r n r n+l 
i = Ct ■ ' 

- F m-i F m a Fn +2 F m , i ?m L "t 1 F n ri+2 a r+Zf™" 

n n— 1 n n— 1 n+1 n n+1 " 



After some calculation, we have 

„, rpm—1 rpm+1 7 77m 77m _. rpm+2 77m 7 T?m+l rpm+1 

Oii'n-l i'n+l ~ a^n+l-^n-1 _ at n+1 1 n-1 ~ a* n-1 *n+l 



77m- 1 7777-1+ 1 77m rpm+2 

1 n 1 n 1 n 1 n 



(85) 



After decoupling, we obtain 



aiC-i KfT " L Jn+iK-i = SF^F^ 1 , (86) 



which is nothing but Eq. (14ip , Thus we conclude that this special case gives a 
discrete potential BSQ equation which has the same r-function to Eq. (|40j) . 
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4. Conclusion 

We have proposed an alternate form of discrete potential BSQ equation. We have 
constructed the bilinear equations and multisoliton solutions for the discrete po- 
tential BSQ equation. The bilinear equations and multisoliton solutions have been 
constructed by one of 3-reductions of the Hirota-Miwa equation. Using the discrete 
potential BSQ equation, we have presented the ultradiscrete potential BSQ equa- 
tion. We have also studied the lattice potential Boussinesq equation of Nijhoff et al. 
using the singularity confinement test. Although the lattice potential Boussinesq 
equation of Nijhoff et al. is in very complicated form, we can find bilinear equations 
easily by using singularity confinement test. We have investigated the relationships 
among our alternate discrete potential BSQ equation, the discrete potential BSQ 
equation of Date et al. and the lattice BSQ equation by Nijhoff et al. 

An interesting problem is to present explicit forms of soliton solutions of the 
ultradiscrete potential BSQ equation. Since 3-reduction condition is very compli- 
cated, it is not easy to see which solutions can survive in the ultradiscrete limit. 
We will address this problem in the near future. 
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